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Abstract 

We consider the dynamics of a harmonic crystal in the half-space with zero boundary 
condition. It is assumed that the initial date is a random function with zero mean, finite 
mean energy density which also satisfies a mixing condition of Rosenblatt or Ibragimov 
type. We study the distribution m of the solution at time t E R . The main result is 
the convergence of /xt to a Gaussian measure as t — > oo which is time stationary with 
a covariance inherited from the initial (in general, non-Gaussian) measure. 
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1 Introduction 



The paper concerns the problems of the long-time convergence to the equilibrium distribution 
for the discrete systems. For one-dimensional chains of harmonic oscillators the problem is 
analyzed in [1, 16]: in [16] - for initial measures which have distinct temperatures to the left 
and to the right, and in [1] - for a more general class of initial measures characterized by a 
mixing condition of Rosenblatt- or Ibragimov- type and which are asymptotically translation- 
invariant to the left and to the right. For many- dimensional harmonic crystals the convergence 
has been proved in [13] for initial measures which are absolutely continuous with respect to 
the canonical Gaussian measure. In [4]-[7] we have started the convergence analysis for 
partial differential equations of hyperbolic type in R d , d > 1 . In [8]-[10] we extended the 
results to harmonic crystals. In the harmonic approximation the crystal is characterized 
by the displacements u(z, t) G R ra , z G Z d , of the crystal atoms from their equilibrium 
positions. The field u(z, t) is governed by a discrete wave equation. In the papers mentioned 
above the lattice dynamics has been studied in the whole space Z d . 

In the present work the dynamics of the harmonic crystals is studied in the half-space 
Z%, d>l, 

u( z ,t) = - Yl ( v ( z ~ z ') -V(z-z'))u(z',t), zeZ d + , teR, (l.i) 

where z := (—z 1 , z) , z = (z 2 , ■ ■ ■ , z d ) G Z^ 1 , with zero boundary condition, 

u(z,t)\ Zl=0 = 0, (1.2) 

and with the initial data 

u(z, 0) = uq(z), u(z, 0) =u 1 (z), z G Z d + . (1.3) 

Here Z+ = {z G Z d : z\ > 0} , V(z) is the interaction (or force) matrix, (Vki(z)) , 
k,l = l,...,n, u(z,t) = (ui(z,t), . . . ,u n (z,t)) , u (z) = (u 01 (z), . . . ,u 0n (z)) G R™ and 
correspondingly for Ui(z) . To coordinate the boundary and initial conditions we suppose 
that Uq(z) = Ui(z) = for Z\ — . 

Denote Y(t) = (Y°(t), Y x (t)) = (u(-,t),u(;t)) , Y = (Y °X) = M0> «i(0) ■ Then 
(1.1)— (1.3) takes the form of the evolution equation 

Y(t)=A+Y(t), teR,zeZ d + , Y°(t)\ zi=0 = 0, Y(0)=Y o . (1.4) 

Here A+ = I ?, I) , with V + u(z) := E (V(z - z') - V(z - z'))u{z') . 

u / Z 'ez d + 

It is assumed that the initial state Y is given by a random element of the Hilbert space 

T~t a ,+ of real sequences, see Definition 2.1 below. The distribution of Yq is a probability 

measure /xo satisfying conditions S1-S4 below. In particular, the initial correlation function 

Q (z,z') is asymptotically translation-invariant as z 1 ,z[ — > +oo (see Condition S2) and the 

measure /x has some mixing properties (see Condition S4). Given t G R, denote by /i t 

the probability measure on H a ,+ giving the distribution of the random solution Y(t) to the 

problem (1.4). 
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Our main result gives the weak convergence of measures [i t on the space 7i a ,+ , with 
a < — d/2 , to a limit measure /Jqc , 

— 7 fXoo as oo, (1.5) 
where /ioo is an equilibrium Gaussian measure on Ti a ,+ ■ This means the convergence 

J f{Y)fH(dY) - 1 /(y) Aioo (dy), t - oo, 

for any bounded continuous functional / on T~t a ,+ ■ 

Explicit formulas for the correlation functions of the measure /xq are given in (3.3)-(3.5). 

The paper is organized as follows. The conditions on the interaction matrix V and the 
initial measure /iq are given in Section 2. The main result is stated in Section 3. Examples of 
harmonic crystals and the initial measures satisfying all conditions imposed are constructed 
in Section 4. The convergence of correlation functions of fit is established in Section 5, the 
compactness of /i t , t > , and the convergence of characteristic functionals of \i t are proved 
in Sections 6 and 7, respectively. 



2 Conditions on the system and the initial measure 

2.1 Dynamics 

Let us assume that 

V(z) = V(z), where z:=(-z u z), z = (z 2 , ■ • . , z d ) G Z d ~\ (2.1) 

Then the solution to the problem (1.1)— (1.3) can be represented as the restriction of the 
solution to the Cauchy problem with the odd initial date on the half-space. More precisely, 
consider the following Cauchy problem for the harmonic crystal in the whole space Z d : 

v(z,t) = - E V(z- z')v(z',t), zeZ d , teR, 

z'eW (2.2) 
v(z,0) = v Q (z), v(z,0) = vi{z), zeZ d . 

Denote X(t) = (X°(t), X\t)) = «, t), v{-, t)) , X = (X°,X*) = (v (-), Vl (-)) . Then (2.2) 
has a form 

X(t) = AX(t), t E R, X(0) = X . (2.3) 

Here A = ^ ^ ) ' W ^ 6re ^ ^ S a convom ^ on operator with the matrix kernel V . 

Let us assume that the initial date X (z) be an odd function w.r.t. Z\ G Z 1 , i.e., 
X (z) = — Xq{z) . Then the solution v(z,t) of (2.2) is also an odd function w.r.t. z\ G Z 1 . 
Let us restrict the solution v(z,t) on the domain Z d + and put u(z,t) = v(z, t)\ Zl >o . Then 
u(z,t) is the solution to the problem (1.4) with the initial date Y (z) = X (z)\ Zl >o ■ 

Assume that the initial date Y of the problem (1.4) belongs to the phase space Tt a ,+ , 
a G R , defined below. 
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Definition 2.1 H a ,+ is the Hilbert space of R n x R n -valued functions of z E endowed 
with the norm 

\\y\\l + = E \Y(z)\ 2 (i + \z\ 2 r<oo. 

zez d + 

In addition it is assumed that the initial date Y = if z 1 — . 
We impose the following conditions E1-E6 on the matrix V . 

El. There exist positive constants C,7 such that ||V(2)|| < Ce^ 2 ! for z e Z d , \\V(z)\\ 
denoting the matrix norm. 

Let V(0) be the Fourier transform of V(z) , with the convention 

V(0) = V(z)e lz -\6e T d , 
zez d 

where " • " stands for the scalar product in Euclidean space R d and T d denotes the d -torus 
R7(27rZ) d . 

E2. V is real and symmetric, i.e., Vi k (—z) — V k i(z)e'R, k,l = \,...,n, z G Z d . 

Both conditions imply that V{6) is a real-analytic Hermitian matrix-valued function in 
6 E T d . 

E3. The matrix V(9) is non- negative definite for every 9 G T d . 
Let us define the Hermitian non-negative definite matrix, 

Q(9) = (V(9)) 1/2 > 0. (2.4) 

£1(0) has the eigenvalues < u\(0) < uo 2 (0) . . . < co s (0) , s < n, and the corresponding 
spectral projections H a (0) with multiplicity r a = trII CT (0). \— > u) a (0) is the a-th band 
function. There are special points in T d , where the bands cross, which means that s and 
r a jump to some other value. Away from such crossing points s and r a are independent of 
. More precisely one has the following lemma. 

Lemma 2.2 (see [8, Lemma 2.2]). Let the conditions El and E2 hold. Then there exists a 
closed subset C* C T d such that we have the following: 

(i) the Lebesgue measure of is zero. 

(ii) For any point 9 G T d \C* there exists a neighborhood 0(Q) such that each band function 
ui a (0) can be chosen as the real- analytic function in 0(0) . 

(Hi) The eigenvalue u a (0) has constant multiplicity in T d \ C* . 
(iv) The spectral decomposition holds, 

n(6) = ^^(6)1^(6), 0eT d \C«, (2.5) 

a=l 

where H a (0) is the orthogonal projection in R™ . U a is a real- analytic function on T d \C* . 
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For 9 G T d \ C* , we denote by Hess (uj a ) the matrix of second partial derivatives. The 
next condition on V is the following: 

E4. Let D a (9) = det (Hess(w ff (0))) . Then D a {9) does not vanish identically on T d \ C* , 

<T=1,...,S. 

Let us denote 

C = {9 G T d : det V{9) = 0} and C CT = {9 G T d \ C, : ,0,(0) =0}, <7 = 1, . . . , s. (2.6) 

Then the Lebesgue measure of C a vanishes, a = 0, 1, s (see [8, Lemma 2.3]). 
The last conditions on V are the following: 

E5. For each a ^ a' , the identities ui a {6) ± uv(#) = const± , 9 G T d \ C* , do not hold with 
const ± 7^ . 

This condition holds trivially in the case n — 1 . 
E6. HV^ 1 ^)!! G L 1 (T d ) . 

If Co = , then ||V^~ 1 (^)|| is bounded and E6 holds trivially. 

Denote by H a the Hilbert space of R™ x R n -valued functions of z G Z d endowed with 
the norm 

ra= e TOi 2 (i+N 2 r<oo. 

zez d 

Proposition 2.3 (see /#, Proposition 2.5]). Let conditions El and E2 /ioZc? ; and choose 
some a G R. Tnen (%) /or any X G H a , there exists a unique solution X(t) G C(H,H a ) 
to the Cauchy problem (2.3). 

(ii) The operator U(t) : X \— > «s continuous in H a ■ 

Corollary 2.4 Let conditions (2.1), El and E2 hold. Then (i) for any Y G T~L a ,+ , there 
exists a unique solution Y(t) G C(R,7i Qi+ ) to ine mixed problem (1-4)- 
(ii) The operator U + (t) : Y \— > F(t) zs continuous in Tl a;+ . 

Proof. Corollary 2.4 follows from Proposition 2.3. Indeed, the solution X(z,t) of (2.3) 
admits the representation 

X(z,t)= £ Qt(z-z')X (z'), (2.7) 
z'ez d 

where the Green function Gt{ z ) has the Fourier representation 

Q t {z) : = F^Jexp (A(9)t)} = (2n)- d J e^ 6 exp (A(9)t) d9 (2.8) 

with 

Therefore, the solution to the problem (1.4) has a form 

y(z,o= £ S t , + MW), ^zf- (2.io) 

where z') := C/t(z - 2') - Q t (z - z!) . Corollary 2.4 follows. ■ 
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2.2 Random initial data and statistical conditions 

Denote by fiQ a Borel probability measure on 7Y Qi + giving the distribution of Y . Expecta- 
tion with respect to /i is denoted by E . 

Assume that the initial measure /x has the following properties. 

51. Y (z) has zero expectation value, E (Y (z)) = 0, z G Z d + . 

For a,b,c G C™ , denote by a <g> b the linear operator (a<S>b)c — a Y7j=i bjCj . 

52. The correlation matrices of the measure /x have a form 

g^^) = ^(^)(8)^V)) = 9^i^i^-^)> z,z'EZ d + , ij = 0,1. (2.11) 

where 

(i) qo(zi, z[, z) — for ^ = or ^ = , 

(ii) lim y ^ +00 q l o(z 1 +y, y, z) = c${z), z = (zi,z) G Z d . Here c$(z) are correlation functions 
of some translation invariant measure v with zero mean value in 7i a . 

Definition 2.5 A measure v is called translation invariant if viT^B) = v{B) , B G B(TL a ) , 
h G Z d , where T h X{z) = X(z - h) , z G Z d . 

53. The measure /x has a finite variance and finite mean energy density, 

e (z) = E (\Y (z)\ 2 + \Y 1 (z)\ 2 )=tr[Q (z,z)+Q 1 1 (z,z)} < e < oo, z e Z d + . (2.12) 

Finally, it is assumed that the measure /x satisfies a mixing condition. To formulate this 
condition, let us denote by a (A) , A C Z d + , the a -algebra in 7i a ,+ generated by Y (z) with 
z E A. Define the Ibragimov mixing coefficient of the probability measure fiQ on 7i at+ by 
the rule (cf. [12, Definition 17.2.2]) 

\fi Q (AnB)-fi (A)^(B)\ 
(p(r) = sup sup — . (2.13) 

ABczj a e <r(A), b e o-(B) /M-dJ 

dist(„4, S) > r Mo(-B) > 

Definition 2.6 A measure fio satisfies the strong uniform Ibragimov mixing condition if 
ip(r) — > as r — > oo . 



oo 



S4. The measure /Iq satisfies the strong uniform Ibragimov mixing condition with 



r d - x ^ j2 {r)dr < oo . (2.14) 

o 

This condition can be considerably weakened (see Remarks 3.4 (i), (ii)). 

3 Main results 

Definition 3.1 (i) We define ji t as the Borel probability measure on H a ,+ which gives the 
distribution of the random solution Y(t) , 

f H(B) = f jL (U+(-t)B), B G B(Ti a ,+), teH. 



5 



(ii) The correlation functions of the measure // t are defined by 

Qi j (z,z') = E^Y^t) ®Yi(z',t)), i,j = 0,1, z,z'eZ d + . (3.1) 

Here Y l (z,t) are the components of the random solution Y{t) = (Y°(-, t), Y 1 ^, t)) to the 
problem (1-4) ■ 

The main result of the paper is the following theorem. 

Theorem A Let d,n > 1 ; a < —d/2, and assume that the conditions (2.1), E1-E6 and 
S1-S4 hold. Then 

(i) the convergence in (1.5) holds. 

(ii) The limit measure //oo is a Gaussian measure on TL a ,+ ■ 

(Hi) The correlation matrices of the measures // t converge to a limit, for i,j — 0, 1 , 

Q?(z, z') = J (Y\z) ® Y\z')) n t {dY) -> Ql{z, z'), f -> oo, z, z' e Z d + . (3.2) 
The correlation matrix Q l ^ {z,z') = (Q^z, z'))}j =0 of the limit measure //oo has a form 

Qoo(z, z') = q^z - z') - qoo(z - z') - q^z - z') + q^z - z'), z,z eZ d + . (3.3) 
Here qoo(z) = q^(z) + q^o(z) , where in the Fourier transform we have 

C(0) = \t, U ^)(^(O) + C(e)q o (0)C(0r)U a (9), (3.4) 

4 er=l 

<U#) = ^E^^i^W) 11 -^)^)^^)-^)^)*)^^), 6eT d \C*,(3.5) 

4 er=l 

H a (0) is the spectral projection from Lemma 2.2 (iv) and 

(iv) The measure //oo is time stationary, i.e., [C/ + (t)]*// 00 = //oo , t e R. 

(v) The group U + (t) is mixing with respect to the measure //oo , i- e -, 

lira J f{U + {t)Y)g{Y)^{dY) = J f{Y)^{dY) J g(Y) ^(dY) 

for any f,g e L 2 (H a:+ , Hoc) ■ 

The assertions (i), (ii) of Theorem A can be deduced from Propositions 3.2 and 3.3 below. 

Proposition 3.2 The family of measures {// t , t G R} is weakly compact on the space H a ,+ 

for any a < — d/2 , and the bounds sup £'||C/ + (i)io||a + < 00 hold. 

t>o 
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Set S = [S'(Z^) ® R n ] 2 , where S(Z d h ) stands for the space of rapidly decreasing real 
sequences. Denote (Y,V)+ = (Y°,^°} + + (Y\^ 1 ) + for Y = (Y°, Y l ) G H a , + and # = 
G S, where 

(Y\V) + = £ Y\z)-V{z), * = 0,1. 
zez^ 

Proposition 3.3 For every $6 5, the characteristic junctionals converge to a Gaussian 
one, 

/e^+^^^expj-iQoo^,*)}, t - oo, (3.7) 
where Qoo is the quadratic form defined as 

Proposition 3.2 ensures the existence of the limit measures of the family {fi t , t G R} , while 
Proposition 3.3 provides the uniqueness. They are proved in Sections 6 and 7, respectively. 
The assertion (Hi) of Theorem A is proved in Section 5, item (iv) follows from (1.5) and item 
(v) can be proved using a method of [11]. 

Remarks 3.4 (i) The uniform Rosenblatt mixing condition [15] also suffices, together with 
a higher power > 2 in the bound (2.12): there exists 5 > such that 

e(\Y °{z)\ 2+s + \Y \z)\ 2+s ) <C<oo, zeZ d + . (3.8) 

Then (2.14) requires a modification: 

r+oo 

/ r d - l a p (r)dr < oo, with p = min(5/(2 + 5), 1/2). 
Jo 

Here a(r) is the Rosenblatt mixing coefficient defined as in (2.13) but without fio(B) in the 
denominator: 

a(r) = sup{a Y (A,B) : A,B C Z d + , dist(„4, B) > r}, 

where 

ay(^,i3) = sup{|/i (An J B)-/io(A)/i ( J B)| : A G 5 G 

Under these modifications, the statements of Theorem A and their proofs remain essentially 
unchanged. 

(ii) The uniform Rosenblatt mixing condition also could be weakened. Let K(z, s) = 

d 

n [zi — s, Zi + s] , where s > , z G Z d , stand for the cube in Z d , if r = Z d \ K(z, s + r) . 

i=l 

Let us define the mixing coefficient cty(r) by the rule 

ai{r) = sup [a Y (K(z, s) n , X r n ) : z G Z^_, < s < l} . 

To prove Theorem A it suffices to assume, together with (3.8), that 

C l K 
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with some constants C, k, k' > . See [2, 3] for a more detailed discussion about the different 
mixing conditions. 

(iii) The condition E5 could be considerably weakened. Namely, it suffices to assume the 
following condition: 

E5' If for some a ^ a' one has 0J a {9) ± u a i(9) = const ± with const-t ^ , then 



p^{o)T^{e)u al {e) P Z{e) = ^ 



(3.9) 



Here p^ a , (9) stand for the matrices 

P^(0):=n,(W(0)lM*), ^T d , a,a' = l,...,s, i,j = 0,1, 

qo (0) are Fourier transforms of the correlation functions qo (z) . 

Note that the condition E5' is fulfilled, for instance, if qo(z) is a covariance matrix of a 
Gibbs measure u on H a , ol < —d/2. Formally the Gibbs measure uq is 

l -^J2(\Mz)\ 2 + Vv (z).v (z)) 
u (dv , dvx) = —e z ^ d JJ dv Q (z)dv 1 (z) , 

where Zp is the normalization factor, (3 = T~ l , T > is an absolute temperature. Let 
us introduce the Gibbs measure u as the Gaussian measure with the correlation matrices 
defined by their Fourier transform as 

cC(0)=TV-\9), c&{9)=T{5 kl )l l=l , = $»{0) = Q. 

Then p%, = for a ^ a' and (3.9) holds. 



4 Examples 

Let us give the examples of the equations (1.1) and measures /j,q which satisfy all conditions 
E1-E6 and S1-S4, respectively. 



4.1 Harmonic crystals 

The conditions E1-E6 in particular are fulfilled in the case of the nearest neighbor crystal, 
i.e., when the interaction matrix V(z) = (V k i(z)) k j =1 has a form: 



Vm(z) 
V kk (z) = < 



for k ^ I, 

-7fc for \z\ 



1, 



2-fk + 'mi for 2; = 0, k — 1, . . . , n, 



(4.1) 



for \z\ > 2, 

with 7^ > , rrik > . Then the equation (1.1) becomes 

u k (z, t) = (7feA L - m 2 k )u k (z, t), k = 1, . . . , n. 
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Here Al stands for the discrete Laplace operator on the lattice Z d , 

A L u(z) : = ^2 (u(z + e) -u(z)). 

e,\e\=l 

Then the eigenvalues of V(9) are 



Cj k (6) = yj 271(1 - cosfli) + ... + 2 7d (l - cos0 d ) + m\ , k = l,...,n. 

These eigenvalues have to be labelled as follows 

£>i(0) = ... = a) ri (0), Q ri+1 (9) = ... = u r2 (9), . .., <Zv._ 1+ i(0) = ... = £v.(0) 

U0i(6) = Cu ri (9) < U0 2 {9) = LU r2 (9) < . . . < uo s (9) = uo rs (9). 

Clearly conditions E1-E5 hold with C* = . In the case all m k > the set C is empty 
and condition E6 holds automatically. Otherwise, if m k = for some k , C = {0} . Then 
E6 is equivalent to the condition uu~ 2 (9) e L 1 (T d ) , which holds if d > 3 . Therefore, the 
conditions E1-E6 hold provided either (i) d > 3 , or (ii) d = 1, 2 and all m k > . 
In the case (4.1) formulas (3.4) and (3.5) can be rewritten as follows. Denote 

, s / 1 if k,l e [r CT _i + l,r CT ] 
^ ((7)= otherwise <r = W, 



Then 



where 



..I s 

Qlki = tE^W. *,j = 0,l, k,l = 1, . . . ,n, 

4 (7=1 



<4q°° + q* 1 - i sign(sin ^H^ 1 - q^°) 



hi ' 



4.2 Gaussian measures 

We consider n — 1 and construct Gaussian initial measures //q satisfying S1-S4. Let us 
define u in H a by the correlation functions c$(z — z') which are zero for i ^ j , while for 

i = 0, 1 , 

q«(0) := F z ^eK{z)\ E L\T d ), <$(9) > 0. (4.2) 

Then by the Minlos theorem there exists a unique Borel Gaussian measure u on 7i a , 
a < —d/2 , with the correlation functions qo (z — z') , because 

\X\\lv Q {dX)= ^(l + |^| 2 r(trq° (0)+trqJ 1 (0))=C(a,rf) J \x(<R(9)+<R(9)) d9 < 00. 

zGZ d T d 

The measure z/ satisfies SI and S3. Let us take a function £ e C(Z) such that 

/ \ f 1, for s > a, 
C( s ) — \ n r ^ n wrt h an a > 0. 
w 0, for s < 



Let us introduce X(z) as a random function in probability space (Ti a , ^o) • Define a Borel 
probability measure /io on H a , + as a distribution of the random function Y (z) = ((zi)X(z) , 
z E Z\ . Then correlation functions of /i are 



Qo(z, z) = <£(z - z')C(*i)C(4), i,j = 0, 1, 



where z = (z±, . . . , Zd) , 2' = (^i, • • • , ^) G , and qg are the correlation functions of the 
measure z/ • The measure /x satisfies SI S3. Further, let us provide, in addition to (4.2), 
that 



Then the mixing condition S4 follows with ip(r) = , r > r$ . For instance, (4.2) and (4.3) 
hold if we set c${z) = f(z 1 )f(z 2 )-. ■ ■•f(zd) , where f(z) = N — \z\ for \z\ < N and f(z) = 
for \z\ > N with N := [r /\fd] (the integer part). Then f(9) = (1 - cos N 8)/(l - cos9) , 
9 E T 1 , and (4.2) holds. 

5 Convergence of correlation functions 
5.1 Bounds for initial covariance 

Definition 5.1 By l p = l p (Z d ) <g> R™ (P + = l p (Z d + ) <g> R n ) , p > 1 , ti>l,we denote tee 
space of sequences f(z) = (fi(z),...,f n (z)) endowed with norm \\f\\ip = ( \f( z )\ p ) ^ 



The next Proposition reflects the mixing property of initial correlation functions. 

Proposition 5.2 (i) Let conditions S1-S4 hold. Then for i,j — 0, 1 , tee following bounds 
hold 



<iS(z)=0, \z\>r . 



(4.3) 





E \Qo( z , z ')\ < C<oo forallzEZ d + 



(5.1) 



Yl \Qo(z,z')\ < C<oo forallz'eZ d + 



(5.2) 



zezi 



Here the constant C does not depend on z, z' G Z| . 
(ii) ^ eC(T d ), i,j = 0,1. 

Proof, (i) By [12, Lemma 17.2.3], conditions SI, S3 and S4 imply 



Ql\z,z')\<Ce ^ 2 (\z-z'\), z,z'eZ d + 



(5.3) 



Hence, condition (2.14) implies (5.1), 



£ \Qi j (z,z')\<Ce Q J2^(\z\)<oc. 



(5.4) 
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(ii) The bound (5.3) and condition S2 imply the following bound: 

\<fi(z)\ <CW /2 (M), zeZ d . (5.5) 

Hence, from (2.14) it follows that qtf(z) E I 1 , what implies qtf E C(T d ) . ■ 

Corollary 5.3 Proposition 5.2 (i) implies, by the Shur lemma, that for any $, ^ E 1+ the 
following bound holds: 

\(Q (z,z>),^z)®y(z>)) + \<c\\n P+ \\nii- 

5.2 Proof of the convergence (3.2) 

From condition (2.1), formulas (2.8) and (2.9) it follows that Gt( z ) — Gt{z) with z = 
(—z 1 , z 2 , ■ ■ ■ , z d ) . Then, by the explicit representation (2.10), the covariance Qt(z, z') can be 
decomposed into a sum of four terms: 

Q t ( z ,z') = y, GtAz,y)Qo(y,y')GfA z '>y') = Rt(z,z')-R t (z,z')-R t (z,z') + R t (z,z'), 

y,y'eZ* 

where 

R t (z,z')-.= ]T Gt(z-y)Qo(y,y')Gl(z'-y'). 

Therefore, (3.2) follows from the following convergence 

Rt(z, z') — >• q^z - z') as t -> oo, z,zEZ d . (5.6) 
To prove (5.6) let us define 



n ( 7 j\ - I Qo(^ z ') for z ^' e Z d + , 

V*\ z ' z )-\ o otherwise. 



First we split the function Q*(z,z') into the following three matrices 

Q + (z,z') := iq (*-*') > (5-7) 

Q-(z,z') : = iqo^-^sign^i), (5.8) 

Q r (z,z') := Q*(z,z')-Q + (z,z')-Q-(z,z'). (5.9) 
Next introduce the matrices 

R?(z,z')= £ (G t {z-y)Q a {y,y')Gl{z' -y')), z,z'EZ d , t > 0, (5.10) 

y,y'&Z d 

for each a = {+, — , r} , and split R t (z, z') into three terms: R t (z, z') = Rf{z, z')+R~ t ~(z ) z') + 
Rl(z,z') . Then the convergence (5.6) follows from the following lemma. 

Lemma 5.4 (i) lim Rf(z,z') = q^(z — z') , z,z' E Z d , with the matrix g+ defined in 

t— +oo 

(34), 

(ii) lim R t (z, z 1 ) = q^iz — z r ) , z, z 1 E Z d , with the matrix defined in (3.5). 

t — > oo 

(Hi) lim R r t (z,z') = 0, z, z' E Z d . 

t— >oo 

This lemma can be proved using the technique of [9, Proposition 7.1]. To justify the main 
idea of the proof we sketch the proof of Lemma 5.4 (i) in Appendix. 
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6 Compactness of measures family 

Proposition 3.2 follows from the bound (6.1) by the Prokhorov compactness theorem [17, 
Lemma II. 3.1] by a method used in [17, Theorem XII. 5. 2], since the embedding H at+ C 7i./3,+ 
is compact if a > (5 . 

Lemma 6.1 Let conditions SI, S3, S4 hold and a < —d/2. Then the following bounds 
hold 

supE\\U + (t)Y \\l + <oc. (6.1) 

Proof. Definition 2.1 implies 

E\\Y(;t)\\l, + = E (1 + \z\ 2 r(trQ t (z,z) +tr Ql\z,zj) < oo. 

Since a < —d/2 , it remains to prove that 

sup sup ||<5t(z, z')\\ < C < oo. 

The representation (2.10) gives 

QVM = E(Y\z,t)®Y\ z ',t))= y: E Gt^y)Q»M)<3i l A^y') 

= (Q (y,y'),& z (y,t)®& zl (y',t)) + , 

where 

®i(y,t) := (gi° + (z,y),gi\(z, y )) 

= (g?(z-y)-g?(z-y),gi\z-y)-gi\z-y)), < = 0,1. 
Note that the Parseval identity, formula (8.2) and condition E6 imply 

\m-,t)\\l = (2vr)- d / \*\(e,t)\ 2 M <C J (\Q?{6)\ 2 + \Gm\ 2 ) d6<C < oo. 

Then Corollary 5.3 gives 

\Qi j (z,z')\ = \(Q (y,y'),<S>i(y,t)®<S>i,(y',t)) + \ < C\\Q),(;t)\\e + WA'Mil < Ci < oo, 
where the constant C\ does not depend on z, z' G , t G R . ■ 



7 Convergence of characteristic functionals 

We derive (3.7) by using the explicit representation (2.10) of the solution Y(t) , the Bernstein 
'room - corridor' technique and a method of [4]-[9]. The method gives a representation of 
(Y(t), as a sum of weakly dependent random variables (see formula (7.6) below). Then 
(3.7) follows from the central limit theorem under a Lindeberg-type condition. The similar 
technique of the proof is applied in [9, Sections 9, 10]. Then we remark only the main steps 
of the proof. 



12 



7.1 Asymptotics of U' + (t)V 

At first, let us evalute of scalar product (Y(t), \&) + . Let us introduce a function ^*(z) as 

r if Zl > o, 

= { 0, if 2i = 0, 
{ -V(z), if ^1 < 0. 

Therefore 

(Y(z, t), *(z))+ = (Y(z, t), *.(z)) + = (Y (z'), Q>(zf, t)) + , (7.1) 

where 

:= i/;(t)*.oo = E &T+(*>* , ) 1 M*) = E 

= (2n)- d j e- iz '- e g; (6)^,(6) d6. (7.2) 

I'd 

Definition 7.1 (%) T/ie critical set C := C U C* U ( Uf C CT ) ('see E4j. 
(iij 5° := G 5 = [S(Z d ) ® R n ] 2 : £(0) = in a neighborhood of C} . 

Note that mes C = (see [9, lemma 7.3]) and it suffices to prove (3.7) for ^ G 5° only. 
For the function &(z,t) the following lemma holds. 

Lemma 7.2 (^c/ Lemma 9.1 from [9]) Let conditions E1-E4 and E6 hold. Then for any 
fixed G 5° ? £/ie following bounds hold: 

(i) su VzeZd Mz,t)\<C r d / 2 . 

(ii) For any p > £/iere eaxsi C p ,7 > such that 

\<S>(z,t)\<C p (l + \z\ + \t\)- p , \z\>-yt. (7.3) 

This lemma follows from (7.2), (8.2), Definition 7.1 and the standard stationary phase 
method. 



7.2 Bernstein's argument 

Let us introduce a 'room - corridor' partition of the half-ball {z G : |z| < 7t} , with 7 
from (7.3). For t > we choose A t and p t G N . Let us choose < 5 < 1 and 

Pt-t 1 - 5 , A t ~-^, t^oo. (7.4) 
log* 

Let us set h t = A t + p< and 

a' = j/i t , V = a* + A t , j — 0,1,2, ... , N t — [(jt)/ht]. 

We call the slabs R{ = {z G Z% : |z| < iVA, a J ' < < &*} the 'rooms', C\ = {z G Z% : 
M < JVt/it, 6* < 21 < a-? +1 } the 'corridors' and L t = {z G Z% : |*| > A^} the 'tail'. Here 



13 



z = (zi, . . . , Zd) , A t is the width of a room, and p t of a corridor. Let us denote by xl the 
indicator of the room R\ , that of the corridor C{ , and rj t that of the tail L t . Then 

N t -1 

where the sum J2t stands for j2 ■ Hence, we get the following Bernstein's type represen- 

3=0 

tat ion: 



<Xo, H; t)) + = E t [( y o, xiH; t)) + + (lo, eM; t)) + ] + (lo, t)) + . 
Let us define the random variables r\ , , l t by 

W = (F ,^$(-,t)) + , c?' = (*o, Z t = (Y , »/*$(.,*))+. 

Therefore, from (7.1) and (7.5) it follows that 

(Y(t), *> + = (y , *(-, *)>+ = EtW + c?) + k- 

Lemma 7.3 Let S1-S4 hold and ty* E S° . The following bounds hold for t > 1 : 

E\rl\ 2 < C(*) A t /t, Vj, 



(7.5) 



(7.6) 



E\ 



< C(*) pt/t, Vj, 



E\l t \ 2 < C p (*)t-r, \/p>0. 

The proof is based on Lemma 7.2 and Proposition 5.2 (i) (see [9, Lemma 9.2]). 

Further, to prove (3.7) we use a version of the central limit theorem developed by Ibrag- 
imov and Linnik. If Q^^, \I/) = 0, the convergence (3.7) follows from (3.2). Thus, we may 
assume that for a given $ * G 5° , 

Qoo(*,*)^0. (7.7) 

At first, we obtain 



|£exp{i<y ,$(-,*)>+}-£oo(*)l = 



+ o(l), £ — > oo. 



This fact follows from Lemma 7.3, convergence (3.2), condition S4 and (7.4) (cf [9, p. 1073- 
1075]). 

Secondly, by the mixing condition S4, we derive that 



Nt-l 



Eexp{iJ2/t}- II Eexp{iri} 



< CN t (p(pt) -> 0, t -> oo. 



Hence, it remains to check that 



JVt-l 

E exp jirf' | — exp 



„J|2 



0, t — > oo. 
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According to the standard statement of the central limit theorem (see, e.g. [14, Theorem 
4.7]), it suffices to verify the Lindeberg condition: V5 > , 



Here <r t = Et-^l r t | 2 > an d ^/ = E{X e f) , where X e is the indicator of the event |/| > e 2 . 
Note that (3.2) and (7.7) imply that a t — > Q 00 (\I / ,\I / ) 7^ , £ — > 00 . Hence it remains to 
verify that 

^2 E £ \rl | 2 — ► 0, t — > 00, for any £ > 0. 
This condition is checked using the technique from [9, section 10]. ■ 

8 Appendix. Outline of the proof of Lemma 5.4 (i) 

Obviously, the assertion of Lemma 5.4 (i) is equivalent to the next proposition. 
Proposition 8.1 Let conditions E1-E6 and S1-S4 hold. Then for any f £ 5, 

Urn {R+(z, z'), (*) ® OO) = (<& (* - A ® *(^))- (8-1) 

Proof. It suffices to prove (8.1) for \1> G S° only. It can be proved similarly as in [9, Lemma 
7.6]. 

At first, let us apply the Fourier transform to the matrix Rf(z,z') defined by (5.10): 

Et(9,9') :=F z ^ e Rt(z,z>) = g t (0)Q + (0,0 , )g?(9'), where Q+(0,0') ■= F z ^ e Q + (z,z>). 
z' ->e' z' -> e' 

From (5.7) it follows that Q + (9, 9') = 5(9 + 9') (2n) d q o (0)/2 . Hence, 

Rj(9,9') = (2^6(9 + 9%mo(0)<3l(-0)- 
Secondly, Gt(9) has a form 

r (o\ I cosfit sinfitfi-M 

WJ-^-ginnfn cosfit J' (8 - 2) 

where fi = Q(9) is the Hermitian matrix defined by (2.4). Let C(9) be defined by (3.6) and 
/ be the identity matrix. Then 

Cj t (9) = cos Sit I + sin QtC (9). (8.3) 

Moreover, by condition E2, Qj(-9) = Q*(9) = cos fit / + sin fit C(9)* . Therefore, 

(Rt(z,z>)Mz)®*(z>)) = (2n)- 2d (Rt(9,9>)M9)®ir(9i)) 

15 



Further, let us choose certain smooth branches of the functions H a (9) and u a {9) to apply 
the stationary phase arguments which require a smoothness in 9 . We choose a finite partition 
of unity 



M 



£> m (0) = l, tfGSUpp^, 



1.5) 



m=l 



where g m are nonnegative functions from C^(T d ) and vanish in a neighborhood of the set 
C defined in Definition 7.1, (i). Further, using (8.5) we rewrite the RHS of (8.4). Applying 
formula (8.3) for Qt{9) , one obtains 

^(v'),$(#^))~E E / ^(^(n.m^wn^w,^^)®^))^, 

l J TO o-,CT' = l Td 

where Rt,aa'(9) stands for the 2n x 2n matrix, 

Rt,aa'(0) = W{^ s M e ) ± ^AO))t[^o(o)i : C(e)c lo (e)C(ey] 



± 



+ 



sm(u; a (e) ±u a ,{9))t [C(0)q o (0) ±qo(0)C(0)*]}. 



(8.6) 



If a = a' , then 

Rt,<ra{9) = 



qo(0) + C(0)q o (0)C(0)*j + - cos (2*^(0)*) [q o (0) - C(0)q o (0)CW 
+^ sin (2w ff (0)f) [C(0)qo(0) + q ( W)*] . (8.7) 



By Lemma 2.2 and the compactness arguments, we choose the eigenvalues u a (9) and the 
matrix H a (9) as real-analytic functions inside the supp g m for every m: we do not mark 
the functions by the index m to not overburden the notations. Let us analyze the Fourier 
integrals with g m . 

At first, note that the identities u a (9)+ui a '(9) = const + or uj a {9)—uj c i{9) = const_ with 
the const ± ^ are impossible by condition E5. Furthermore, the oscillatory integrals with 
uj a {9) ±u a r(9) ^ const vanish as t — > oo . Hence, only the integrals with u a (9) — uj a i{9) = 
contribute to the limit, since uj a {9) +cu a i(9) = would imply u a {9) = ou a >(9) = which is 
impossible by E4. By formulas (8.6) and (8.7), one obtains 

(Rt(z,z')^(z)®**{z')) 

= (2n)- d J2i2i I ^)(n^)[qo^) + c^)qo(0)c^r]n CT ^) + ...,^)®¥^))^ 

(7 = 1 4 



= (2vr)- d + 



where " ..." stands for the oscillatory integrals which contain cos(cj ct (#) ± uj a >(9))t and 
sin(a; o -(0) ±w (7 /(0))t with w cr (6 l ) ±u; (7 /(6') ^ const. The oscillatory integrals converge to zero by 
the Lebesgue-Riemann theorem since all the integrands in ' ... ' are summable, and V(u a (9)± 
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oj a i{9)) = only on the set of the Lebesgue measure zero. The summability follows from 
Proposition 5.2 (ii) and E6 (if C 7^ 0) since the matrices H a (9) are bounded. The zero 
measure follows similarly to Lemma 2.2 (i) since u a {9) ± u a >{9) ^ const . Lemma 5.4 (i) is 
proved. ■ 
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